Hydrodynamics of superfluids confined in blocked rings and wedges 
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Motivated by many recent experimental studies of non-classical rotational inertia (NCRI) in 
superfluid and supersolid samples, we present a study of the hydrodynamics of a superfluid confined 
in the two-dimensional region (equivalent to a long cylinder) between two concentric arcs of radii b 
and a {b < a) subtending an angle /3, with < /3 < 2tt. The case (5 = 2-n corresponds to a blocked 
ring. We discuss the methodology to compute the NCRI effects, and calculate these effects both for 
small angular velocities, when no vortices are present, and in the presence of a vortex. We find that, 
for a blocked ring, the NCRI effect is small, and that therefore there will be a large discontinuity 
in the moment of inertia associated with blocking or unblocking circular paths. For blocked wedges 
(6 = 0) with /3 > TT, we find an unexpected divergence of the velocity at the origin, which implies the 
presence of either a region of normal fluid or a vortex for any nonzero value of the angular velocity. 
Implications of our results for experiments on "supersolid" behavior in solid *He are discussed. A 
number of mathematical issues are pointed out and resolved. 

PACS numbers: 47.37.-|-q, 47.32.Ef 
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I. INTRODUCTION 

Flow without dissipation is the defining feature of su- 
perfluidity. Because of this property the moment of in- 
ertia of a vessel containing a superfluid is different from 
(smaller than) that when the liquid is in the normal state. 
This effect is largest in the absence of vortices, when su- 
perfluid flow is irrotational. The difference between the 
moments of inertia when the liquid, confined by bound- 
ary conditions, is in the normal and superfluid states is 
known as the "non-classical rotational inertia" (NCRI). 
The occurrence of NCRI is often used as an experimental 
signature of superfluidity. Superfluid hydrodynamics and 
the resulting NCRI have been studied extensively [ij in 
the past for simple geometries, such as spherical, cylindri- 
cal or rectangular containers rotating about a symmetry 
axis. Because of several recent developments, some of 
which are briefly discussed below, it has become neces- 
sary to understand the properties of flow of superfluids in 
enclosures of more complicated geometry. These provide 
the motivation for our present study. 

Recent observations 0, i, i, i, 1 S 1 of NCRI in 
torsional oscillation experiments on solid '*He have been 
interpreted as the occurrence of a "supersohd" phase. 
This interpretation of the experimental results is con- 
troversial. There is experimental and theoreti- 
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cal [lOl evidence suggesting that the observed NCRI 
is due to superfluidity along crystalline defects such as 
grain boundaries in a polycrystalline sample and net- 
works of dislocation lines. Since these extended defects 
form complex disordered structures, calculations of the 
flow properties and the rotational inertia of a superfluid 
confined in irregular-shaped channels are necessary for 
a quantitative assessment of whether this mechanism is 
the correct explanation of the observed results. In this 
context, it is important to examine whether the super- 
fluid component can flow along continuous closed paths 
in the sample. Since the geometry of the network of 
defects would depend on thermodynamic variables such 
as temperature and pressure, and on the cell geometry, 
the availability of such paths would also depend on these 
parameters and conditions. Thus, an understanding of 
the dependence of the NCRI on such variables requires, 
for example, a calculation of how the NCRI arising from 
a blocked ring of superfluid changes as the blockage is 
removed. To check whether the observed NCRI is due 
to the occurrence of extended superfluidity, the NCRI 
of samples in which the solid ^He is confined in the an- 
nular region between two concentric cylinders has been 
measured 0, 0] in the presence of a barrier in the annu- 
lus that prevents possible flow of the superfluid along a 
closed path surrounding the rotation axis (the common 
axis of the cylinders). The NCRI observed under these 
conditions is found to be much smaller than that for sam- 
ples in which the artificial block is not present. The cal- 
culation just mentioned is obviously relevant for a quan- 
titative understanding of the results of such experiments. 
Finally, an understanding of experimental results [1, 0, Q 
on the dependence of the NCRI on the frequency of tor- 
sional oscillations requires a theoretical analysis of vortex 
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formation and critical velocity in superfluids confined in 
irregular-shaped channels. 

Our study is also partly motivated by the recent ex- 
plosion of activity in experimental and theoretical stud- 
ies of superfluidity and other quanturn phenomena in 
trapped, ultracold atomic systems [H, Also, there 
have been many experimental studies of the flow prop- 
erties and NCRI of superfluids confined in porous me- 
dia such as yyc or g lass and containers packed with fine 
powder [3, Ea, floT [l7l| . The first experimental observa- 
tion of "supersolid" behavior was in a torsional os- 
cillator experiment on solid ^He confined in vycor glass. 
Since the pores in these systems have complex geometry, 
it is necessary to work out the hydrodynamics of super- 
fluids in irregular-shaped channels in order to understand 
the results of these experiments in quantitative detail. 

Thus, we study here the hydrodynamics of a super- 
fluid confined in a two-dimensional region between two 
concentric circular arcs, each of which subtends an angle 
(3 at their common center. The annular region between 
the two arcs is bounded on two sides by straight walls 
along the radial direction. Thus, the special case with 
/3 = 27r corresponds to a ring that is blocked by a wall 
placed perpendicular to its inner and outer peripheries. 
This two-dimensional geometry corresponds, neglecting 
edge effects, to that used in many experiments on super- 
solid behavior in ^He where the helium is confined in the 
annular region between two concentric cylinders, under 
the assumption that the cylinders are long enough and 
the confined system is homogeneous along the cylinder 
axis. In the limit of vanishing inner radius, the geometry 
we study corresponds to that of a wedge with opening 
angle (3. The limit /3 = 27r in this case represents a cir- 
cular container with a straight blocking wall extending 
from the center of the circle to its periphery. 

We assume throughout the paper that the fluid is in- 
compressible, which is appropriate for superfluid Helium. 
We flrst consider the case where there are no vortices (so 
that the superfluid flow is irrotational), and solve the hy- 
drodynamic equation for the velocity fleld for rotation 
about an axis perpendicular to the plane of the system 
and passing through the common center of the arcs that 
form its boundary. The sample geometry is reflected in 
the boundary conditions for the velocity fleld. For in- 
compressible and irrotational flow, the velocity fleld can 
be expressed in terms of either a scalar or a vector poten- 
tial (stream function), analogous to those in electromag- 
netic theory, both of which satisfy the Laplace equation 
with appropriate boundary conditions. The scalar poten- 
tial method is simpler, and leads to series that converge 
rapidly. We have used this method to obtain the veloc- 
ity field for (3 = 2-K and (3 — n. For a general value of 
/3, however, the stream function method, although more 
difficult in that it leads to series that are not convergent, 
but Borel summable, is more powerful. We have there- 
fore used it to obtain the velocity field for arbitrary (3. 
We present analytic results for the velocity field and the 
moment of inertia for arbitrary values of the inner and 



outer radii and the opening angle (3. We also derive a 
simple "parallel axis" theorem that relates the moment 
of inertia for rotation about any axis perpendicular to the 
plane of the system to the calculated value for rotation 
about an axis passing through the center of mass. 

In the context of experimental observations of NCRI 
in solid '*Hc, the most important result of our study is 
about the NCRI of a blocked ring. When the ring is 
blocked, the superfluid can not flow through it. How- 
ever, due to the irrotational nature of superfluid flow, 
the moment of inertia is smaller than that for rigid-body 
rotation. Therefore, the drop in the moment of inertia 
when the block is removed (the superfluid does not con- 
tribute to the moment of inertia when there is no block) 
is less than the rigid-body value. Our calculations show 
that the moment of inertia of a blocked ring whose width 
is small compared to its radius is very close to its moment 
of inertia for rigid rotation, so that unblocking the ring 
(i.e. the opening up of a closed path) produces a large 
drop in the moment of inertia (nearly equal to its rigid- 
rotation value), which would show up in an experiment 
as a relatively large value of the NCRI. Thus, the onset of 
NCRI in experiments on solid *He may correspond to the 
unblocking of large closed paths in the network of defects 
along which the superfluid component is supposed to ex- 
ist. Our results for the NCRI of a superfluid conflned 
in a blocked ring can be compared directly with those of 
experiments 0, |3| in which the NCRI of solid ^He con- 
flned in an annular cell is measured both in the presence 
and in the absence of a barrier that blocks flow around 
the annular channel. We show that our results, when 
combined with accurate measurements of the NCRI, can 
provide valuable information about the structure of the 
superfluid network in solid "^He, and discuss the validity 
of our hydrodynamic description for superfluid flow in 
narrow channels such as those along crystalline defects 
in solid ''He. 

An interesting new result of our calculation is that the 
velocity fleld for a wedge with f3 > ir diverges at the tip 
of the wedge for any nonzero value of the angular velocity 
fi. This means that the implicit assumption that the ve- 
locity field nowhere exceeds the Landau critical velocity 
is in principle mathematically incorrect for these wedges: 
for any nonzero value of il, there must be a region near 
the tip where the liquid is in the normal state. We show 
that the size of the region where this occurs is too small 
to have any measurable consequence in ''He experiments 
performed with usual geometries. This divergence of the 
velocity can be removed by the presence of a single vor- 
tex. We calculate the position of this vortex and the 
rotational inertia in its presence. Our calculations un- 
cover also several interesting mathematical issues and we 
indicate ways of addressing them. Some of these were 
also present in earlier studies [ij of superfluid hydrody- 
namics, while some are new. We discuss these questions 
as they appear throughout the paper. 

Whether vortices appear or not is in general deter- 
mined by the free-energy cost of creating a vortex. We 
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will show that for typical experiments on ^He, vortices do 
not occur for sufficiently small angular velocities. How- 
ever, as pointed out in Ref. [ij, states with vortices 
present will have, at sufficiently larger values of the angu- 
lar velocity, a lower free energy than the vortex-free state. 
We calculate the critical angular velocity for vortex nu- 
cleation which turns out, for typical *He samples, to be in 
the experimentally important range of angular velocities. 
We show how the rotational inertia is modified by these 
vortex excitations. These results are relevant for under- 
standing the experimentally observed dependence of the 
NCRI of "supersolid" '^He on the frequency of torsional 
oscihations [1, Sill- 

The rest of this paper is organized as follows. In sec- 
tionini we describe in detail our calculations. We present 
first two alternative methods of calculating the velocity 
field in the vortex free case, and discuss the results ob- 
tained for this field and the moment of inertia. We com- 
pare our results for the NCRI with those of experiments 
on solid *He in blocked annular geometry and point out 
other implications of our results for experimental studies 
of superfluidity. We then explain how to include vortices 
in our description and calculate the critical angular ve- 
locity for vortex nucleation. A summary of our results is 
presented in the concluding section ITlTl 



II. RESULTS 
A. Formulation of the problem 

We consider, as explained above, superffuid ffow in an 
ideal cylinder, long enough in the z direction so that 
edge effects are negligible and the problem quasi two- 
dimensional. The cross sections of the cylinders that we 
will consider will be bounded by two concentric circular 
arcs of radii a and h (with a > b) and encompassing an 
angle that we will call (3. In the limit b — the shape 
of this cross section is that of a circular wedge. We will 
consider all values of /3, < < 2tt. It must be empha- 
sized that the case f3 — 2tt is not the same as that of a 
ring, since a boundary along a radius still exists. 

In the absence of vortices (the generalization to the 
case when vortices are present will be discussed below) 
the superfluid velocity field v(r) for an incompressible 
fiuid satisfies: 

V-v(r)=0 (2.1a) 
V X v(r) = 0. (2.1b) 

The boundary condition corresponding to superffuid ro- 
tation around some center O with uniform angular ve- 
locity H is that [H the normal component of the fluid's 
velocity at the boundary must equal the normal compo- 
nent of the rigid-body velocity x r at that point. That 
is, the component of v(r) along the outward normal n to 
the boundary must equal, at any point on the boundary. 



the component of Jl x r along n at that point: 

v(r) • n = (r2 X r) • n (2.2) 

were r is a vector from O to a point on the boundary. 
The point O is not necessarily the center of mass of the 
system: in general we will take it to be, for reasons of 
obvious computational convenience, the center of the arc 
or arcs that are part of the boundaries of our system. 

There are two obvious ways to solve Eqs. (I2.ip . The 
first is to introduce a scalar potential V{r) such that 
v(r) ~ W{r). In that case T^(r) satisfies the Laplace 
equation, V^V^(r) = 0, and Eq. (|2.2p is a Neumann 
boundary condition on V. Alternatively, one can intro- 
duce a stream function ^'(r) such that: 

v., = -d^/dy (2.3a) 
Vy = d'ii/dx, (2.3b) 

where one can think of ^ as the z component of a vector 
potential [v(r) = —V x (z\I'(r))]. It is obvious that ^'(r) 
also satisfies the Laplace equation, V^5'(r) — 0. Now, 
however, the boundary conditions are of the Dirichlet 
form [l[: at any point in the boundary, 

*(r) = ^Qr^. (2.4) 

It turns out, as we will see, that for certain special val- 
ues of P such as TT and 27r, the scalar potential method 
is much simpler to use and leads to expressions for v(r) 
in the form of rapidly convergent series which are very 
convenient. However, for other values of /?, this method 
becomes rather awkward. The stream function method 
on the other hand can be used for any value of /?, but the 
resulting expressions involve asymptotic series. These 
are, however, Borel summable and agree with the re- 
sults obtained from V(r) in the cases where the scalar 
potential method works well. For this reason, we will 
first present here results obtained from ^(r) for (3 = 2tt 
and (3 — 7: and then consider the general case using the 
stream function. 

Once the velocity field is obtained, the angular momen- 
tum (and hence the moment of inertia) can be calculated 
by straightforward integration of the velocity field. In 
this way, the depletion of the moment of inertia from its 
rigid body value is obtained. In general our origin O is 
not the center of mass (COM) of the system: therefore 
it is important to discuss an interesting property of the 
nature of the parallel axis theorem shift in the superfluid 
case. If one considers the moment of inertia of the super- 
fluid with respect to the COM, Igp'^^ one finds, of course, 
that it is always smaller than that of the corresponding 
rigid object (RO) of the same shape and density, I^q^^ ■ 
Indeed, for the case of a circle Ig'p^ vanishes. With re- 
spect to an arbitrary origin O one has for the superfluid 
a total moment of inertia Igp = Igp^^ + Igp where the 
last term is the parallel axis shift. The key point here 
is that this shift is the same as that for the rigid object. 
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One has: 



jPA _ J 



PA 
RO- 



(2.5) 



The proof of this theorem is very simple: the problem, 
as defined by the above equations and boundary condi- 
tions, is linear. If one shifts the origin from the COM to 
a point a distance R away from it, the velocity field of 
the boundaries shifts to v = (r + R) x $1. In view of this, 
the linearity of the problem, and the boundary condition 
Eq. (|2.2[) , the solution of the shifted problem is the veloc- 
ity field computed with respect to rotations around the 
COM, plus a uniform velocity field R x SI. This second 
field trivially satisfies the equations and takes care of the 
additional term in the boundary condition. But it is triv- 
ial to verify that such a constant field leads simply to a 
parallel axis theorem shift in the moment of inertia equal 
to that for the corresponding rigid object. This applies 
irrespective of the shape of the object: it is not limited to 
the wedge shapes considered here. It is straightforward 
to check by direct calculation that it applies, for exam- 
ple, to the ellipsoidal shapes of Ref. [l|. This theorem has 
physical consequences: since the parallel axis shift cannot 
be "depleted" from its RO value by the superfiuid fiow, 
in general the fractional depletion of Isf will always be 
largest when the rotation is around the COM. 



radius). For this case, one can very simply use the scalar 
potential method. We write, in polar coordinates: 

V{r,(t>) = a™r"/2 gin(TO(/)/2)-H^ fo^r™/^ cos(to0/2). 

m>l m>l 

(2.6) 

With the radial wall set along the = direction, the 
azimuthal component of the velocity. 



(^'<^)= E ya™r-™/2-' cos(m0/2) 



m>l 



(2-7) 

- E ^&™r'"/2-isin(m0/2) 



must equal fir at ^ = 0. This immediately tells us that 
all the am vanish except 04, which equals 11/2. The radial 
component is then: 

Vr{r,(t>) = flrsin(20) + ^ ^6„r"/2-i cos(m0/2) 

m>l 

(2.8) 

At r = a we have Vr = and from this one obtains that 
all the hn with even n are zero while for odd n: 



B. Scalar potential method for P — 2tv and (3 = n 



32r2o 



Trn{n'^ — 16)a"/^ 



(2.9) 



To illustrate the results, let us first turn to the simplest 
case where (3 = 27r, 6 = (a circle with a wall along its 



From these and Eqs. (|2.7p and (|2.8p we have the final 
result for the velocity field: 



J 



16i7a ■r-^ 
Vr[r, (p) — Qr sin(20) H > 

77- / J 



n>0, n odd 



n2 - 16 



v^{r,(j)) — r2rcos(2(/i) 



E 



n>0, n odd 

I 



— 16 



cos{n(f>/2) 



sm{nip/2) 



(2.10a) 
(2.10b) 



where p = r/a. 

Two remarks are needed about these simple results: 
first, the series involved are very rapidly convergent. Sec- 
ond, the velocity components have a square root sin- 
gularity at the origin. Mathematically, the singularity 
is integrable, and allows for the formal calculation of 
the moment of inertia. Physically, the relevant num- 
ber is the value of r at which the velocity would ex- 
ceed the Landau critical velocity Vc- For liquid ''He, 
Vc ~ 2.5 X IC'cm/s [1^, and in typical experiments on 
supersolid behavior, the maximum value of f2 is less than 
0.1s~' (see for example, 0, 0|). This would mean than 
only at values of r/a around 10~^' would Vc be exceeded. 
Such small values of r would not have any experimentally 



measurable consequence (the hydrodynamic description 
we use would not even apply to such length scales). Also 
this divergence is not present for nonzero values of the in- 
ner radius &, and the inner radius is finite (of order lO"'^ 
cm) in torsion and rotation experiments. Thus, this di- 
vergence is not important for ^He. This divergence may 
have observable consequences in Bose-Einstein conden- 
sates (BEC) in cold atomic systems [H, [l^, although 
our incompressibility and uniform density assumptions 
are not applicable to BEC in cold atomic systems, where 
the high compressibility and the confining potential in- 
troduces substantial variations in the density. We show 
later that the divergence discussed above is present in 
blocked wedges for all values of /3 greater than tt. The 
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effects of tliis divergence are discussed in sections III Dl 
andUm 

The angular momentum is obtained by integration of 
rv^ over the sample and the moment of inertia is just the 
ratio of the angular momentum and the angular velocity 
fl. We will use units in which the areal mass density is 
unity. We obtain the result: 



128a4 



SF 



y 

^ - 16)(n + 4) 

n>0, n odd ^ ' 



(2.11) 



which, after numerically evaluating the rapidly conver- 
gent series, gives Isf = 0.693a'*. Thus we have for this 
obstructed circle: 



Isf 
Iro 



0.441. 



(2.12) 



The same method can be used at /3 = tt. In that case 
the only significant difference is that in the expression for 
V{r) one must write: 



V{r, 0) 



a„ir"^ sm{m<j>) + ^ bmr"^ cos(m0). 

m > 1 ?n > 1 

(2.13) 

As before, all the coefhcients a„ are determined from the 
boundary conditions on at (/> = and (j) = n. Both 
are satisfied if all a„ vanish except ai = f2/2. The 6„ 
are determined then from the boundary condition on Vr- 
The result for the velocity field is: 



.(r, (/)) = nr sin(2(/)) + > 



P" 7cos(n0) 



TT — ' 71^—4 
n>0, n odd 



Vfi,{r, (p) — fir cos(2(j)) - ^ p" -sm{n(j)). 



(2.14a) 
(2.14b) 



n>0, n odd 

r 



The series are again convergent, and now the previously 
found integrable singularity at the origin is absent. The 
moment of inertia with respect to the origin is: 



16a4 



n>0, 71 odd 



1 



4)(n- 



(2.15) 



Numerically, we have Isf = 0.488a'* which gives a ratio 
Isf/Iro = 0.621, a value higher than that for the circle. 
However, we must recall that in this case O is not the 
COM and that (as shown above) there is no reduction in 
the parallel axis term so that from the point of view of 
the COM the reduction must be larger. Indeed one finds 
that: 



tCOM 
^SF 
rCOM 
^RO 



= 0.41, 



(2.16) 



which is actually a little less than that for the circle. 

One can see that it is awkward to extend this simple 
procedure to other values of (3. If one sets for example 



P — tt/2 and doubles again the angles and powers in 
the expression for V{r) one finds that it is not possible 
to satisfy the boundary condition for at (j) = and 
= 7r/2 from a single term in the first sum (the a„ co- 
efficients) in the potential. Similar difficulties are found 
at e.g. (3 = 37r/2. Although these difficulties should not 
be unsurmountable, we will instead use the stream func- 
tion method in the general case and deal appropriately 
there with the mathematical difficulties associated with 
the asymptotic series that then result. 

However, one can easily generalize this simple proce- 
dure, for the above values of /3, to the physically more 
relevant case where b > 0. We will consider here the im- 
portant case of an obstructed ring, P = 27r. In that case 
one simply has to add to the potential in Eq. (|2.6p the 
appropriate negative powers of r. The coefficients are 
then found from the boundary conditions on at r = a 
and r = b. One then obtains the velocity fields: 



J 



, , . , , 1651a X ^ / , 1 
Vr{r,(j)) ^ ilapain{2(j)) > cos(n0/2)- — - 

TT ^-^ (1 — C")(l 



n>0, n odd 



(1 -C")(7l2 - 16) 



W0(r, (j)) — flap cos 20 



E sin(n</./2)- ^ 



n>0, n odd 



(l-c")(n2-16) [' 



(2.17a) 



(2.17b) 
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FIG. 1: The velocity field for a blocked ring with c = 0.5. The 
first panel shows the relative strengths of the velocity field as 
a function of position. The second panel is the radial compo- 
nent (in units of fla) plotted vs p = r/a at azimuthal angles 
(from bottom to top) 7r/16, 7r/8, 7r/4, 77r/4, ISvr/S, 3l7r/16. 
The third panel, in the same units, shows the azimuthal com- 
ponent of the velocity vs. at p = 0.6, 0.75, 0.9. 



where c = 6/a < 1, /„(c) = 1 - c"/^+^ and .g„(c) = 
c" - c"/2+2. Plots of the fields given by Eqs. (|TT7|) are 
shown in Fig. [T] All the plots in the figure are for c = 
0.5, a value in the region where, as we shall see below, 
NCRI effects are found to be largest. In the first panel, 
the vector field is displayed in two dimensions over the 
entire sample. The units of velocity are arbitrary, but the 
overall pattern of the field is then clearly shown. In the 
second and third panels we show a plot of Vr (in units of 



ria) vs r (in units of a) at several values of the azimuthal 
angle (p and a plot, in the same units, of Vtj, vs (j) at several 
values of r. One can see that the boundary conditions 
are satisfied. 
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FIG. 2: Moment of inertia of an obstructed ring in terms of 
its aspect ratio c = b/a. In the top panel the ratio R of Isf 
(Eq. (|2.18p l to the rigid body value is plotted, while in the 
bottom panel we plot Isf itself, in units such that a = 1. The 
maxima in the two plots are at different values of c. 



The moment of inertia of the superfluid blocked ring 



is: 



128a4 



iSF 



TT ^ ?i(n2 - 16)(1 - c") 

n>0, n odd ^ ' 



E 

1 

n + 4' 



n - 4 



(2.18) 



The behavior of this quantity as a function of aspect ratio 
c is well worth noting. In the first panel of Fig. [5] we plot 
the ratio R = Isf/Iro for a blocked ring of aspect ratio 
c, vs. c. As noted above, the value for c = (blocked cir- 
cle) would, strictly speaking, have to be corrected, but 
the range of c affected by this is negligible. The ratio 
R increases very quickly with c: at c = 1/2 it already 
reaches 0.875 while at c = 0.75 it exceeds 97%. We see, 



therefore, that a narrow superfluid circular channel rotat- 
ing about its center behaves essentially like a rigid body 
when it is blocked. Since, when unblocked, its moment of 
inertia vanishes, we see that in such a channel there will 
be a sharp discontinuity in / as it is blocked or unblocked. 
In a sample containing a number of such channels, dis- 
continuities or glitches in / will occur as the channels are 
blocked or unblocked. As c — > 1, i? — > 1 and the un- 
blocking would drop R from one to zero, the maximum 
amount. One should recall, however, that / vanishes at 
c — 1 for both the superfluid and the rigid body. In an 
experimental situation one would measure the difference 
in / with the channel blocked and unblocked which is 
Is F itself. This quantity has a broad maximum centered 
around c « 0.52 as one can see in the second panel of 
Fig. [21 There we plot Isf itself in units such that a is 
unity. From this plot one can see that the important 
experimental contribution would come from a range of 
rings with c values in the region 0.2 through 0.8. 



C. Stream function method for arbitrary {3 



As discussed in section III Al the velocity field can be 
written in terms of a stream function ^'(r) that satisfies 
the Laplace equation with Dirichlet boundary conditions 
(see Eqs. (ESm^])). Following Ref. fi*], the general solu- 
tion for ^'(r) for arbitrary (3 can be written as 



*(r) = j dl'r'^n' ■ V'G(r',r), 



(2.19) 



where the line integral / dl' is over the boundary of the 
system, n' is a unit vector along the outward normal to 
the boundary, and G(r, r') is the Green's function for the 
Laplacian operator, satisfying the equation 



V2G(r,r') =(5(r-r'). 



(2.20) 



and the boundary conditions G(r, r') = for all r on 
the boundary of the system. Thus, \E'(r) and hence, the 
velocity field, can be obtained from Eq. (|2.19p once an 
expression for the Green's function, satisfying Eqs. (|2.20|) 
and its boundary condition is obtained. 

As in section Hi B1 we first consider, for simplicity, the 
case 6 = 0, which corresponds to a wedge of radius a 
and opening angle (3. The Green's function in this case 
is easily obtained [20| to be 



G{r, <t>; r', ^') - -- E ^ ^ Mnncf>/P) Mnncb' / ft), (2.21) 

n=l \^> / 

where r> (r<) is the larger (smaller) one of the two radial coordinates r and r' . Using this in Eq. (I2.19p . we obtain 
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the following expression for the stream function ^'(r): 



1 /rx'^^Z/S 
n \a 



(2.22) 



n>0, n odd 

The radial and azimuthal components of the velocity field, obtained from (f)) through Egs. 12.31 are given by 



Tir ^ — ' \ 3 

n>0, n odd 



7/3 -2 



Va/ a2/ n\a' 



(2.23a) 



'^0('', 0) 



2rJa2 ^ 

TT ^ 

n>0, n odd 



sin(n7r0//3) x 



2r n7r2//32 
a? n27r2//32 - 4 



Va/ ' -"2^2 /fl2 



/3r 



n27r2//32 



(2.23b) 



Calculation of the velocity field for /? = 7r/2 requires some 
care because the denominators of some of the terms in 
Eqs. (P?^ go to zero for /3 = 7r/2 and n = 1. The nu- 
merators also vanish for these values of /3 and n, so that 
finite contributions that vary smoothly with /3 across 7r/2 
are obtained for the velocity components. Similar behav- 
ior is found for /3 = 37r/2 for which the n — term in 
the denominators in Eqs. (|2.23p vanishes. These results 
also exhibit, for /? > tt, a singularity 
which can be readily seen from Eqs. (|2.23p to arise from 
the n — \ term in the sum. This is in agreement with 
what we found from the scalar potential method. As 
discussed in detail in the previous subsection, this diver- 
gence is not physically relevant for '*He, but may have 
observable consequences in experiments on cold atomic 
systems. Its possible physical effects are discussed in sec- 
tions |TTD] and [nEl This singularity is always integrable: 
therefore, the angular momentum of the superfluid about 
the origin (tip of the wedge) is easily calculated for all (i 
using these expressions for the velocity components. The 
result for the moment of inertia about O is 

v-^ 1 (n-K \ 1 
'--V E -(- + 4)^;^^^^. (2.24) 

n>0, n odd ^ ^ ' ^ ' ' ' 

For the case (3 — 2ti, the moment of inertia about the 
origin is given by the infinite series 

Is.{P-2^) = - E -7^- (2-25) 

n>0, n odd ^ ' 

This infinite series appears to be different from the one in 
Eq. (|2.1ip which was obtained using the scalar potential 
method. In particular, the series in Eq. (|2.25|) converges 
more slowly than the one in Eq. (|2.1ip . However, it can 
easily be shown that these two expressions for the mo- 
ment of inertia are mathematically identical. We have 
also checked that a similar situation applies when the re- 
sults for the moment of inertia obtained from Eqs. (|2.23p 
for P = TT are compared to those obtained in the preced- 
ing section using the scalar potential method. 



However, the situation is much more complicated 
when, instead of comparing the moments of inertia, one 
compares directly the velocity fields obtained by the two 
methods. In this case it is not sufficient to add or sub- 
tract a series that converges to zero. The reason is that 
while the series in Eqs. (|2.10p converge for all angles (j) 
and for any r ^ 0, those in Eqs. ((2?23)) and ((2?22)) do 
not. This question is related to other technical difficul- 
ties with the result (|2.22p , and in general with the stream 
function method, which we will further address below. 

The moment of inertia of the wedge for rigid-body ro- 
tation about O is Ino — /3fl*/4, and its moment of in- 
ertia for rigid-body rotation about its COM is given by 
Ir8^' - Iro - with = 8a4 sin2(/3/2)/(9/3). Us- 
ing these results and Eq. (|2.24[) . we have calculated the 
ratios Isf/Iro and Isf^^ /^ro'^^ functions of the an- 
gle p. The results are shown in Fig. [31 These ratios are of 
course less than unity, the level of suppression being given 
by the NCRI effect. In the figure we see that this frac- 
tional suppression is always larger in the COM frame, 
that is, Isf/Iro is always higher than Isf^' I^ro^' ^ 
cept of course aX P — 2ti where the two are the same. 
This is in agreement with the theorem proved at the end 
of Sec. lira It is interesting that the ratio Isf^' I^ro^' 
is not a monotonic function of /? - it exhibits a minimum 
ai P = it/2. 

A representative plot of the velocity field for a wedge 
with P — (7/8)27r is shown in Fig. [H The velocity vector 
field is plotted in arbitrary relative units, as in the first 
panel of Fig. [T] It is instructive to compare that panel 
with Fig. m In the earlier case we have c = 0.5 whereas 
in Fig. m we have a wedge, c = 0. The rise in the absolute 
value of the velocity as r — > can now be seen. On the 
other hand, the behavior of of function of 6 is 

clearly very similar: it follows from the second panel of 
Fig. [H that Vr is very small except for angles near the 
radial boundaries, and this is clearly the case also for 
this c — Q wedge. The behavior of with (p is also quite 
similar. 

We now return to the technical difficulties with the 
general solution for the velocity field obtained above via 
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FIG. 3: The ratios Isf/Iro (upper curve), and Isf'^'/Iro^' 
(lower curve) for a superfluid wedge as a function of the open- 
ing angle /3, < /3 < 2-k. Isf is calculated from Eq. (|2.24p . 



the stream function. As noted in section III Al the quan- 
tity \l/(r, (j>) should be equal to flr'^ /2 at all points on the 
boundary, and the physical velocity field should satisfy 
the boundary conditions Vtf,{r, (p) = rO. for = 0, /3 and 
Vr{r, (f) — Q for r — a. It is easily seen from Eqs. (I2.22[) 
and (I2.23bp that both ^'(r, (/)) and Vff,{r,(f) vanish for 
(/i = and (j) — P (since sm{mT(p / /3) — for these values 
of 4>) . Thus the boundary condition on the radii appears 
to be violated even though the construction of the vector 
potential via the Green's function would seem to ensure 
that it will not be. As to Eq. (|2.23ap for the radial com- 
ponent of the velocity, it can be written as 



-5S 
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FIG. 4: Plots of the velocity field inside the wedge for (3 — 
(7/8)27r. This should be compared with the first panel of 

Fig.m 



Vr{r,(l)) = 



Ecosimrcf)/ (3)- 
1 



n>0, n odd 



(D 



J. ^ mr/P 

a/ 



E 

n>0, n odd 



cos(n7r(/)//3). (2.26) 



While the first term on the right-hand side of this equa- 
tion vanishes for r = a, the second term does not. Thus, 
this component also appears not to satisfy the required 
boundary conditions. Numerically, however, we have 
found that these quantities do approach values consistent 
with the required boundary conditions as the boundaries 
are approached from inside, but there is a discontinuity 
as the boundary is approached and the values exactly at 
the boundaries do not satisfy the boundary conditions. 
This does not affect the calculated values of the angu- 



lar momentum and the moment of inertia because these 
quantities are not sensitive to the values of the velocity 
components exactly at the boundary. 

However, this numerical argument is not fully satisfac- 
tory. Fortunately there are better ones. First, one can 
see that this behavior is associated with the nonconver- 
gence of the series. The last term in Eq. (|2.26p . for ex- 
ample, is not merely nonzero: the series that it contains 
is not convergent while that in the first term is. Indeed 
the rearrangement of terms leading from Eq. (|2.23ap to 
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Eq. (|2.26p isolates just this nonconvergent part. However, 
by rewriting the cosines in terms of exponentials one can 
verify that the series in the last term of Eq. (|2.26p is 
Borel summablel^?] (and also Euler summable) with the 
result being zero. With this proviso, Eq. (|2.26[) satisfies 
the boundary condition analytically. Similar arguments 
can be made for ^Pji and for the azimuthal component of 
the velocity. 

This mathematical problem can also be solved by re- 
defining the stream function as 



where again the first term in the right side is the previous 
result, in this case Eq (j2.23bp . The added quantity is 
zero at all points inside the wedge, and is equal to Qr for 
(j) = 0, P, so that the required boundary conditions for 
these values of <p are now satisfied. The equation for Vr 
becomes 



Vr{r, <p) Vr{r, <p) + —- 

n>0, n odd 



cos(n7r<?!)//3). (2.29) 



X! ~ sin(n7r0//3) - ^ 

n>0, n odd 



(2.27) 

where the first term in the right side is that given by 
Eq. (|2.22p . The second term in the right side, which 
is subtracted from the old expression, is zero for all 
points inside the wedge [2lj, and is equal to — r2r^/2 for 
(j) = 0, p. Therefore, the subtraction of this quantity does 
not affect the behavior of ^'(r, (j)) inside the wedge (where 
it still satisfies the Laplace equation). At the same time, 
the redefined ^'(r, 0) satisfies the required boundary con- 
dition for ip — 0, (3. The new term leads to the following 
additional terms in and Vr'. 



E ■ 

n>0, n odd 



■ sin(n7r0//3) 



4 

(2.28) 



The new term, added to Eq. (|2.23ap . cancels the "offend- 
ing" second term in Eq. (|2.26[) . so that the re-defined Vr 
satisfies the required boundary condition at r = a. 

A similar problem with boundary conditions is also 
present in the solution given in Ref. Ij for the velocity 
field inside a cylinder with a rectangular cross section. 
The expression for the stream function given in Eq. (62) 
of Ref. [Ij does not in fact satisfy the required boundary 
conditions posed there at all points on the boundary. As 
in the case considered here, this does not affect the results 
for the calculated physical quantities in Ref. [l| , and this 
mathematical problem can be cured by the addition of a 
term similar to the one considered above. 

The above calculations can be modified readily to treat 
a superfluid confined in the annular region between two 
concentric arcs with radii a and b {a > b). The Green's 
function in this case has the form 



G{r,<l>;r',c^') = -"E -^^(^y^p^ [rT^' ~ ^] x - ^ j sin(n.0//3) sin(n.07/?)- (2-30) 



In this case one does not have to worry about the behavior as r — > 0. Asymptotic series in the summations over n 
are again encountered and handled as in the preceding case. Using this in Eq. (|2.19p . the stream function ^'(r, 0), 
and from it, the radial and tangential components of the velocity are obtained. We skip the long expressions for these 
quantities and quote the final result for the moment of inertia about the origin: 



16a4 ^r^ 1 



^ n>0, n odd " ' 



+ 4 ,^ 2x„ 1 ^ _4^^_, ^ 22;„ 



(1 - c^) - 5-{(l + c%\ + c^-O - 4c^c-} 



(2.31) 



Here, a;„ = nn/P, c — b/a, and Iro — /3(a^ — ^"')/4 is 
the moment of inertia for rigid-body rotation. We have 
checked that this expression reduces to that in Eq. (|2.24p 
for 6 = 0, and to that in Eq. (|TT5)) for /3 = 27r. In 
Fig. [51 we show results for the NCRI in an annular wedge, 
as obtained from Eq. (|2.3ip . The plots are the same 
as in Fig. [3] except that now we have c — 0.5, in other 
words, the fields are as in Fig. [1] Again, the fractional 
suppression is larger, as it must be, in the COM and it 
exhibits a maximum as a function of 



The results derived above have a direct relevance to 
torsional oscillator experiments on solid ''He 0, 0] in 
which the helium is confined in the annular region be- 
tween two concentric cylinders and the NCRI is measured 
both in the presence and in the absence of a barrier that 
prevents flow around the annulus. If the NCRI in the 
absence of the barrier is due to superflow along a closed 
channel surrounding the rotation axis (the common axis 
of the inner and outer cylinders), then the measured value 
of the NCRI when the barrier is not present should be 
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(lower curve) for an annular wedge (Eq. (12.311) ) plotted as a 
function of the opening angle /3, < /3 < 27r, at a fixed value 
of c = 0.5 



(A/)opon — PsIbo where ps is the supersolid fraction and 
Iro the rigid-body moment of inertia of the channel of 
flow about the rotation axis. The NCRI in the presence of 
the barrier should be given by (A/)ciosed = Ps {Irq-Isf) 
where, if this channel is approximately circular, Isf is 
the moment of inertia of a blocked superfluid ring calcu- 
lated above. Thus, the ratio R' = (A/)cioscd/(A/)opon 
should be equal to {Iro — Isf)/Iro = 1 — i?, where 
R = Isf /Iro depends (see Eq. p38)) and Fig. ^ on 
the value of c = b/a. If the superfluid component were 
distributed homogeneously throughout the sample, then 
a and b would be the outer and inner radii of the an- 
nular cell. Whether this is the case can be determined 
by comparing the experimentally measured value of R' 
with (1 — i?o), where Rq is the value of R obtained from 
Eq. (|2.18p using these values of a and b. If the superfluid 
is instead confined in a channel (or in several separate 
channels) with width substantially smaller than that of 
the annular cell, then R' should be smaller than (1 — Rq) 
because R increases as the width of the ring is decreased 
(see Fig. [2]). 

In the experiment of Ref. [2], a = 0.75 cm and b = 
0.64 cm, so that (1 - Rq) = 0.00817. The experimen- 
tal value of R' is 0.015 which is within a factor of two 
of (1 — Rq) but, surprisingly, it is higher. However, the 
value of (A/)open appropriate for the blocked cell was 
evaluated from the results of a different experiment us- 
ing another cell, so that the quoted value of R' may not 
be very accurate. Also, a value of R' larger than (1 — Rq) 
may be rationalized by assuming that the sample con- 
tains a large number of narrow superfluid channels, most 
of which do not form closed paths around the annulus 
(i.e. have /3 < 2tt). These "naturally blocked" channels 
make small contributions to the net sample NCRI. These 
contributions are not strongly affected by the imposition 
of the external barrier, which can change the value of /3 
for the channels it intersects: our calculation shows that 



R = Isf I Iro for narrow annular wedges with f3 < 2tt is 
rather insensitive to /3. Since these channels contribute 
almost equally to (A/)opon and (A/)ciosod, the value of 
the ratio R' would increase. 

More recently, both (A/)opcn and (A/)ciosod have been 
measured using the same cell [7]. In this experiment, two 
cells, both with a = 0.794 cm and b = 0.787 and 0.745 cm 
were used. In both cases the NCRI in the blocked config- 
uration was found to be smaller than the resolution of the 
experiment. This is consistent with our calculated values 
of (1 - Rq) which are 2.9x10"^ and 1.3x10"^, respec- 
tively. Although the measurements are not sufficiently 
accurate to provide more detailed information about the 
channels of superflow, it is clear that more accurate mea- 
surements of R' for samples with different a and b, com- 
bined with the results of our calculations, would be very 
useful for elucidating the geometry of superfluid channels 
in solid ''He. 

If the superfluid channels are very narrow, the valid- 
ity of the hydrodynamic description used here (and else- 
where 0) might be questioned. However, recent numer- 
ical studies [lO, UH indicate that the diameter of the su- 
perfluid region near the core of a dislocation and the 
width of the superfluid layer along a grain boundary are 
of the order of a few nanometers (^ 10 interparticle spac- 
ings). These values of the superfluid layer width are likely 
to be lower bounds, since superfluid channels of such 
small lateral dimensions can not explain the relatively 
large superfluid density measured in recent torsional os- 
cillation experiments . It has been suggested [1^ that 
the effective lateral dimension of the superfluid region 
near a crystalline defect may be larger due to a kind of 
"proximity effect", as in superconductors. Also, stud- 
ies [23| of the thermodynamics of a system of interact- 
ing vortex lines in type-II superconductors, which can be 
mapped to the zero temperature quantum mechanics of a 
two-dimensional system of interacting bosons, show that 
the width of grain boundaries can exceed 15-20 interpar- 
ticle spacings in some cases. A hydrodynamic descrip- 
tion should be valid if the width of the typical superfluid 
regions is of order ~ 10 interparticle spacings or more: 
this has been well-established quantitatively in several 
numerical studies of the flow pro perties of classical liq- 
uids through narrow channels [25l . [2^ . The same should 
to be true for superfluid ^He because its coherence length 
is very small. 

A related effect that needs to be considered if the su- 
perfluid channel along a crystal defect is very narrow is 
the modulation of the density of the superfluid due to the 
potential arising from the surrounding crystalline region. 
We expect our calculations to be valid in the presence 
of such density modulations because the hydrodynamic 
equation for a rotating superfluid derived (for low angu- 
lar speed) in a recent study [l^l in which superfluidity 
is assumed to coexist with a periodic modulation of the 
density (Eq. (8) of Ref. 27]) is identical to that used in 
our calculation. 
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D. Formation of vortices in a wedge with P > n 

As noted above, the velocity field obtained from a cal- 
culation in which it is assumed to be irrotational exhibits 
a divergence as r ^ for a wedge with (3 > tt. Thus Vc 
must be exceeded near r = 0, implying that either there 
is a region of normal fluid near the tip of the wedge, or 
a vortex is present in the system. As we have indicated, 
this issue is unimportant in the torsional oscillation ex- 
periments because the region of normal fluid near the tip 
would be unobservably small for experimentally relevant 
parameter values. It is, however, interesting to inquire 
about the behavior in the general case. We show here 
that this divergence in the velocity field is eliminated by 



the introduction of a single vortex. 

From symmetry, the vortex must be located along the 
line = /3/2. Let the position of the vortex be (r„, /?/2). 
The presence of a vortex of circulation k{= h/m, where 
h is Planck's constant and m is the mass of a parti- 
cle of the fluid) at {r',(j>') leads to an additional term, 
KG{r,(j3;r',(j)') in the expression for the stream function 
^'(r', (f>) where G(r, (f>; r' , (j)') is the Green's function given 
in Eq. (|2.2ip (see Section 3 of Ref. [l| for a derivation 
of this result). This additional term in ^'(r, 0) (with 
r' = r„, (j)' = f3/2) leads to the following additional term 
in the expression for the radial component of the velocity 
near r = 0: 



oo / n7r//3 \ 

vr{r, c^) = t;°(r, 0) + ^ E '^""^^ " cos(n^</>//3) sin(n^/2) ^ v', + (2.32) 

n=l \rv « 



where v^{r,(j>) is the curl- free result as given by 
Eqs. (|2.29p . The n ~ 1 part of the additional term can- 
cels the divergent n — 1 contribution of the previous 
expression if 



1 



1 



(2.33) 



It is easy to check that the divergence in the expression 
for the azimuthal component of the velocity is also re- 
moved if this condition is satisfied. Defining {ry/ay^^ = 
^, the solution of Eq. f^g)) is ^ = [^/A + r]^ -r]]/2, where. 



K(4-7r2//32) 



> 0. 



(2.34) 



One sees that ^ has the nice property that < ^ < 1 for 
any value of rj. The value of f changes from 1 to as 
the dimensionless parameter 7 = Q,a^/K increases from 
zero to a large value, i.e. the vortex moves inward from 
the rim of the wedge to its tip as the angular velocity 
increases. 

Using the expressions for the radial and tangential 
components of the velocity in the presence of a vortex, 
the total angular momentum of the superfluid can be 
calculated. The presence of the vortex increases the an- 
gular momentum about the origin by the amount Ly and 
the moment of inertia for rotation about the origin by 
Iv — Ly/^- Using the result for the vortex position, 
this can be written as 



640^ 



[4 - 7r2//32][(a/r„)V/3 _ (r„/a)V/3] 



n>0, n odd 



n(4 - n^TT^//^^) 



a 



In the presence of the vortex, the moment of inertia about 
the origin is {Isf + Iv) where Isf is given by Eq. (|2.24p 
and Iv is given by the equation above. The value of 
Ty/a to be used in this equation is given by the solution 
of Eq. (|2.33p . Since the vortex position Vy depends on 
the angular speed fl, the value of ly also depends on il. 

Although the divergence in the velocity field at small r 
is eliminated by the introduction of a vortex, the free en- 
ergy of the state with this vortex is not necessarily lower 
than that of the vortex-free state with a small region 



of normal fluid near ?■ = 0. Specifically, in experimen- 
tal situations (e.g. in experiments on solid ''He discussed 
above) where the dimensions of the region of normal fluid 
are extremely small, the free energy cost of creating the 
normal region is negligible and the free energy cost of 
creating a vortex is the deciding factor in determining 
whether a vortex will be present. We therefore calculate, 
in the following subsection, the free energy of a state with 
a single vortex. 
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FIG. 6: The critical angular velocity for vortex nucleation in 
a ring {j3 = 2tt). Here the critical value of the parameter 7 
(i.e. Qio^ /k) is plotted as a function of c. The circles are 
numerical results, connected by straight dashed lines. The 
increase at larger c shows that the nucleation of vortices is 
unfavorable in that case. 



E. Free energy of a vortex and critical angular 
velocity for vortex nucleation 

In the free energy calculation, we consider the general 
case of a ring with b ^ 0. The angular speed fii at which 



nucleation of a first vortex will occur can be determined 
from free energy considerations. The free energy F is 
given [l| in terms of the energy E and the angular mo- 
mentum L as: 

F = E-m (2.35) 

We will denote here with a subscript the quantities F, 
E and L in the vortex-free state, and with a 1 subscript 
those in the presence of one vortex. As stated in the 
preceding subsection, the stream function in the presence 
of a vortex is: 

«'i(r) = *o(r) + k,G{y, r') = ^-o + l-i (2.36) 

where G'(r,r') is the Green's function given in Eq. ()2.30p 
and r' is the vortex position with coordinates r', 0'. From 
symmetry considerations 4>' = (3/2 and the equilibrium 
radial position of the vortex, r' = ry, is to be determined 
from free energy minimization. The velocity field and the 
angular momentum in the presence of a vortex can be 
readily obtained from the stream function of Eq. (|2.36p . 
The angular momentum is given by 

Li^Lo + na^C (2.37) 

where the dimensionless quantity C has the following ex- 
pression: 



TT ^ — ' n 4: — xi 1 — c^-^" 



with Xn = rnr/p. 

It is not hard to see explicitly that G(r,r') has, as 
expected, a logarithmic singularity at r', so that we can 
write: 

G(r, r') = ^ ln(|r - v'\/a) + g{r, r') (2.38) 

where a is the radius of the vortex core and g(r,r'), 
the nonsingular part of the Green's function, satisfies 



the Laplace equation. As shown in Ref. [l| (see also 
Ref. |28|), the energy in the presence of a vortex can 
be written as 

El = ^L^n + ^KQr'^ - iK*o(r') - ^^^gir', r'). (2.39) 

After some algebra, the nonsingular part of the Green's 
function appearing in Eq. (j2.39[) is obtained as 



5(r',r') ^ ^In (^) - i ^ "r^t^c^" " (^'/«)^^" - W-')^^"], (2-40) 

^ ^ ^ n>0, n odd 

I 



where a;„ — mr/(3. Using Eqs. (|2.37p . (|2.38p . (|2.39p and (r',/3/2) may be obtained. The results depend on the 
p.40p . the free energy in the presence of a vortex at vortex core size, via the logarithmic dependence on a/ a 
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mentioned above. One then minimizes Fi with respect 
to r' to obtain its optimal value and compares -Fi 
and Fq to find the overall equilibrium state. This de- 
pends on the value of f2 and, for sufficiently small fl, it 
is the vortex-free state, while for > fii the one-vortex 
state first becomes favorable. In practice these calcula- 
tions can be done only numerically, but the computations 
are not difficult. The relevant dimensionless parameter is 
the quantity 7 = ^la^/n defined in the preceding section. 
This parameter is the ratio of the characteristic scale, fta, 
of the velocity field v" due to the rotation alone, and the 
scale of the additional velocity field due to the vortex, 
which is k/o. One needs also to input the value of a/ a 
for which we take the physically reasonable value of 10^ 
appropriate for liquid ^He. 

Results for Oi computed for a blocked annular ring 
(/3 — 27t) are given in Fig. [6l There we plot the crit- 
ical value of 7 vs the aspect ratio c. We see that at 
reasonably small or intermediate values of c the critical 
value of 7 is in the range 10-50 corresponding to angular 
speeds in the general range of 10~"'^/s, which is in the 
experimentally relevant region. At large values of c this 
quantity increases, refiecting that the system is behaving 
more like a rigid body, in which case the formation of vor- 
tices is obviously less favorable. A similar trend was seen 
for progressively flatter ellipsoids in Ref. [l|. This implies 
that one need not worry about the formation of vortices 
in narrow blocked rings and wedges while estimating the 
contribution of these objects to the NCRI of the system. 

In Fig. [71 we show the texture of the velocity fleld 
due to the nucleated vortex alone at c = 0.5 and at a 
value of 7 slightly higher than its critical value, which at 
this value of c is 71 ~ 20 (see Fig. [5]). The calculated 
optimal position of the vortex at these values of 7 and 
c is r^/a = 0.74. This position is marked by a (blue) 
circle in the plot. The flelds in this flgure should be 
combined with those in the top panel of Fig. [T] One 
should recall that both plots are in arbitrary units, so 
that before plotting the combined field one should divide 
the fields in Fig. [7] by 7 ~ 20 to take into account their 
overall smaller relative scale. If that were done, however, 
then the plot would be very hard to distinguish with the 
naked eye from that in the top panel of Fig. [T] 

The moment of inertia of a ring in the presence of a 
nucleated vortex may be calculated from Eqs. (|2.37p and 
(|2.38p . The results deviate from those obtained for the 
vortex-free state only by a correction of order I/7. For 
il > r^i this is therefore significant only at small values 
of c. At c ^ we find for example that, at /3 = 2tt, the 
moment of inertia of a blocked wedge (c = 0) increases by 
about 8.3% as a vortex is nucleated at = fii, and the 
increase in the moment of inertia due to the nucleation 
of a vortex becomes less than 1% for c > 0.33. 

The optimal value, r.^, of the radial coordinate of the 
vortex obtained from free-energy minimization is quite 
different from the value for which the velocity due to 
the vortex cancels the mathematical singularity at r = 
found in wedges with /3 > tt. This implies that the veloc- 
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FIG. 7: (Color online) Fields produced by a nucleated vor- 
tex in an obstructed ring with c — 0.5, at = f2i. Only 
the fields produced by the vortex are included. Its position 
(marked by a (blue) open dot) is at the optimal value (see 
text) Vv/a = 0.74. The total flow is the sum of that shown 
in this flgure, weighed by a factor of I/7, and that in the top 
panel of Fig.[T] Because 7 is rather large, the result would be 
hard to distinguish from that shown in Fig. [1] 

ity field would formally diverge at r = in such systems 
even when a vortex is present at the position correspond- 
ing to the minimum of the free energy. As noted above, 
this mathematical singularity does not have any physical 
consequence in usual experiments on ^He. However, this 
interplay between the requirements of keeping the veloc- 
ity below the Landau critical value and minimizing the 
free energy may lead to nontrivial behaviors in other ex- 
perimentally accessible situations such as Bose-Einstein 
condensates in cold atomic systems. 



III. SUMMARY AND DISCUSSION 

We have calculated here the velocity fields of a super- 
fluid sample in a cylindrical wedge, or ring-wedge geom- 
etry. We have used two different methods to solve the 
relevant hydrodynamic equations both in the absence of 
vortices and when vortices are present. From the re- 
sulting velocity fields, we have derived formulas for the 
moment of inertia, and therefore for the NCRI effect in 
these geometries. 

Physically, the most important of our results is that 
the NCRI effect is most prominent for relatively narrow 
rings. Our calculations show that the moment of inertia 
of a blocked narrow ring is very close to the rigid-body 
value unless the width of the ring is a large fraction of 
its outer radius. Since the moment of inertia of a su- 
perfluid ring for rotation about its center is zero when 
it is unblocked (at least for small fi), one should see a 
considerable change in the NCRI when approximately 
circular superfluid channels in a sample are obstructed 
or unobstructed. The fractional change in the moment 
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of inertia as a ring is unblocked (defined relative to the 
moment of inertia of the ring for rigid-body rotation) is 
maximum when the rotation axis passes through the cen- 
ter of the ring. In that case, this ratio approaches unity 
very quickly as the aspect ratio c of the ring is increased 
toward one (see Fig. [21 top panel), and this ratio has a 
value close to 0.44 as c ^ 0. The magnitude of the change 
in the rotational inertia upon blocking/unblocking does 
not depend on the location of the axis of rotation. For a 
fixed value of the outer radius a, the magnitude of this 
change is maximum when the aspect ratio c is close to 
0.52 (see Fig. [5J bottom panel). This maximum is very 
broad. For an annular superfluid wedge, the moment of 
inertia about an axis passing through its tip is close to 
the rigid-body value if the opening angle /? is small, and 
it decreases as /3 is increased (see Figs. |3] and [5]). 

The results summarized above are for the case where 
there are no vortices, so that the velocity field is irrota- 
tional. Since one expects vortices to be nucleated as the 
rotational speed in increased, we have used a free-energy 
criterion to determine the critical angular speed for the 
nucleation of a vortex in the system. We find that in 
standard "supersolid" experiments the relevant range of 
geometries and speeds includes both the parameter re- 
gion where vortices are absent and that where nucleated 
vortices exist. For a fixed value of /3 = 27r (ring geom- 
etry), the critical angular speed increases rapidly as the 
aspect ratio c is increased above about 0.5 (see Fig. [6]). 
Also, the increase in the moment of inertia due to the 
nucleation of a vortex is rather small (less than 10%) in 
all cases. These observations imply that the results men- 



tioned above for a narrow ring without vortices remain, 
for "^He, valid for relatively large values of the angular 
velocity. 

Mathematically, a number of relevant results have been 
uncovered and emphasized. There are a number of tech- 
nical difhculties in the calculation of the velocity fields, 
leading to non-convergent series and singularities. How- 
ever, the singularities are integrable and the series are 
Borel summable, so that there is no difficulty in calculat- 
ing physical quantities such as the angular momentum 
and the kinetic energy. We also point out the occur- 
rence of a mathematical singularity in the velocity field 
in wedges (but not in rings) with (3 > tt and discuss pos- 
sible effects of this divergence. This singularity turns out 
to have no measurable consequence in experimental stud- 
ies of "^He, but may be relevant in studies of cold atomic 
systems confined in wedge-shaped traps. 

In general, the ideas and methods developed here can 
be used in other geometries. We believe that the results 
and techniques presented here can be very useful in un- 
derstanding not only NCRI phenomena in "supersolid" 
helium, but also superflow in confined geometries and in 
finite systems. Work in which we apply these ideas to 
study the NCRI effect in realistic models of grain bound- 
ary networks is in progress. 
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